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Abstract 



A variant of the complex Ginzburg-Landau equation is used to investi- 
gate the frequency locking phenomena in spatially extended systems. With 
appropriate parameter values, a variety of frequency-locked patterns including 
flats, 7r fronts, labyrinths and 2tt/3 fronts emerge. We show that in spatially 
extended systems, frequency locking can be enhanced or suppressed by diffu- 
sive coupling. Novel patterns such as chaotically bursting domains and target 

patterns are also observed during the transition to locking. 
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Frequency locking of nonlinear oscillators is a well known phenomenon, and has been 
studied extensively however, frequency locking in spatially extended systems is not 

well understood. Recently, frequency locking has been observed experimentally in the light 
sensitive Belousov-Zhabotinsky reaction |||4J]. The system is oscillatory, and shows spiral 
patterns without external forcing. When periodically illuminated with light pulses, however, 
it exhibits frequency- locked patterns such as flats, ir fronts, labyrinths and 2tt/3 fronts, 
depending on the forcing frequency. It was also reported that these resonances have the 
structure similar to Arnol'd tongues ||. 

For the theoretical investigation of frequency locking in spatially extended systems, we 
consider the forced complex Ginzburg-Landau equation f6HT3| 

d t ijj = (1 - iv)^ + (-1 + ib)\if;\ 2 ip + (1 + id)V 2 i> + (1) 

where ip(r, t) is a complex order parameter, and u, b, d, c are real. 

A theoretical study on the Eq. ([[]) showed that resonant forcing can exhibit rich patterns 
as varying the external frequency ||. The same equation has been studied for parametric 
resonance cases in the system such as Faraday waves J7|,Q and vertically oscillated granular 
layers ||. However these systems are turbulent || or stationary j7|-[J] without forcing. For 
the oscillatory cases, the studies have been focused on the front dynamics in 2:1 jlO| and 



In : 1 ||nHl3|[ resonance cases. More works are required to understand the frequency locking 
phenomena in spatially extended system. 

In this letter, we analyze the Arnol'd tongue structures of spatiotemporal patterns for 
1:1,2:1 and 3:1 resonance cases. Especially, we describe how diffusive couplings affect the 
frequency locking. 

Firstly, we will explain the model and parameter values. Eq. ([I]) with c = 0, v = 
is a well-known complex Ginzburg-Landau equation, which describes oscillatory media and 



shows stably rotating spiral patterns when 1 — bd > [14-16 



dtifj = ij) + (1 + id)V 2 i) + (-1 + ib)\i)\ 2 i). (2) 
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In this case, the model is considered as the amplitude equation of chemical oscillations 
ho = ipe twt + c.c., where w is the frequency of a spiral rotation in real experiments. 

Introducing a local amplitude A(f, t) and a local phase 4>{r,t) by ip(r,t) = 
A(r,t)exp(i(j)(r,t)), Eq. (fj) reduces to the following phase equation p| , |i5|| : 

8 t 4> = b + (1 - bd) V 2 - (b + d) ( V0) 2 . (3) 

It has a oscillatory solution = (b — (b + d)k 2 )t + k ■ r + <po. Diffusion coefficient 1 — bd 
must be positive for stable phase dynamics. Therefore, two constraints on parameter values, 
1 — bd > and b + d < 0, are required from the experimental facts that without forcing there 
are stably rotating spiral patterns and its frequency is higher than that of a homogeneous 
system 0. Specifically, we choose b = — 0.5, d = —1.4. 

Now, consider when the spirals are forced with external frequency w e = n(w + u), where 
n is an irreducible integer fraction and v is a small deviation from a rational multiple of 
the natural frequency w. In this case, ip is taken as the complex amplitude of chemical 
oscillation h = ipe i{w+u)t + c.c, so that stationary ip may indicate the entrainment of system 
to the external frequency w e . Correspondence between h and h = ipe lwt + c.c. when the 
forcing is absent gives the term —iuip in Eq. ([!]). Considering the invariance of the system 
under a translation t — > t + ^ yields the resonant forcing term cip n 1 in Eq. ([[]) 

Next, we will obtain the Arnol'd tongue structure of Eq. ([J) and study the characteristic 
of resonant patterns. Eq. ([!]) contains two control parameters v (measuring frequency 
mismatch) and c (forcing amplitude). In this parameter space, we seek the stationary 
solutions of Eq. ([]]), which means frequency locking. Without diffusion the linear stability 
analysis of the stationary state yields the Arnol'd tongue of a single oscillator (given as solid 
lines in Fig.[l|). 

There are n equivalent stable stationary solutions for n : 1 resonances. These solutions 
have the same amplitudes, but are phase-shifted with respect to each other. In spatially 
extended systems, different oscillators may lock to n different states to make 27r/n front 
solutions. 
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To obtain the resonance regions of spatiotemporal patterns, we integrate Eq. (JJ) numer- 
ically. System size is 100 by 100 with periodic boundary conditions. The initial condition 
is a set of Gaussian distributed random numbers of mean zero, and variance 0.25. The 
numerical results are summarized in Fig.|l|; the typical patterns are presented in Fig.0. 

When the frequency locking does not occur, basic patterns are spirals. At 1:1 resonance, 
resonant patterns are homogeneous flats (FigJ^(a)). At 2:1 resonance, there are two resonant 
patterns: 7r fronts (Fig.^|(b)) and labyrinthine patterns (Fig.^(d)). Fig.|l|(b) shows that 
labyrinthine patterns arise at higher frequencies than 7r fronts. Under some conditions it 
fronts move until only one domain remains. These front dynamics and the nonequilibrium 
Ising-Bloch transition(NIB) have been studied by previous researchers |5|, [10|JI7| . At 3:1 
resonance, we can see the 27r/3 fronts as in Fig||(c). But the time series at one typical 
point show that these fronts are in fact not stationary, but move slowly(Fig.[3](a)). The 
periodic nature of time series indicates the motion of the fronts. The flat region of time series 
corresponds to the time interval during which each domain passes. Sometimes moving fronts 
merge and eventually collapse into flat patterns. Ref . |]TU[ shows that chiral interfaces move 



in nonequilibrium systems. In this respect, observed 27r/3 front motion can be explained 
from the nonzero chirality of the interfaces (Figj^(b)). 

Next, we discuss how diffusion affects frequency locking and Arnol'd tongue structures 
(See Fig.p. The regions indicated by IF, 2F, 3F and LABY, which corresponds to flat, n 
fronts, 27r/3 fronts and labyrinths, respectively, are frequency locked. The Arnol'd tongue 
structure in Eq. (HD almost coincides with that of a single oscillator, but some deviations 
are found. For example, the tip of Arnol'd tongue is smoothed so that locking begins at 
a nonzero forcing amplitude. It is rather surprising that frequency-locked labyrinths arise 
at parameter values at which a single oscillator without coupling would not be locked. It 
means that diffusion induces frequency locking. One can also see that at lower u, the locking 
is somewhat suppressed. (Note that the regions indicated by BD are not frequency-locked.) 

To understand the enhancement and suppression of locking, we consider the phase dy- 
namics of Eq. ([[]). Substituting ip(r,t) = A(r,t) exp(i(f>(r,t)) into Eq. (U) and assuming 



that the amplitude is slaved to the phase, we get the following phase equation for n = 2: 

d t <j) = (b - u) + (1 - bd)V 2 (j) -(b + d)(V(j)) 2 + cVTTV 2 ' cos(20 + 5), (4) 
where tan(<5) = 1/6. 

Eq. (^|) has an oscillatory solution <p = (b — v — (b + d)k 2 )t + k ■ f + for c = 0. 
The stationary solution for nonzero c means the entrainment of the system to the forcing 
frequency: frequency locking. There are two effects on the frequency of system: forcing and 
nonlinear coupling fl9[]. Forcing, represented by cyl + b 2 cos(20 + S), compensates the term 
(b — v) and entrains the system. Thus, without the nonlinear coupling —{b + d)(V<p) 2 , a 
frequency- locked stationary solution occurs for c > \b — v\j\/\ + b 2 . In other words, forcing 
increases the frequency for b — v < and decreases the frequency for b — v > 0. However, 
nonlinear coupling always increases the frequency of the system by — (b + d)k 2 0]. Hence, 
for high z/, forcing and nonlinear coupling effects are cooperating and frequency locking is 
enhanced. But for low u, two effects are competing and locking is suppressed. 

These effects also explain the following transition routes to locking as increasing the 
external forcing amplitude c. (See Fig.[TJ) 

low v : spiral — > chaotically bursting domain — > frequency-locked front 
high v : spiral — > labyrinth (n = 2), target (n = 3) — > frequency-locked front 
At low z/, the turbulent states, characterized by chaotic nucleation and annihilation of 
domains, arise between spirals and locked states(Fig.|2|(f)). As explained above, the diffusion 
effect against locking increases with (V0) 2 at low v. Near the defect core, (V</>) 2 is large 
and frequency locking is suppressed. However, far from the defect core, (V0) 2 is small and 
domains are frequency-locked. As the defects and the domains move, the dynamics become 
quite complex. The chaotic time series of Re(?/>) at one typical point, shown in Fig.^(a), 
consists of flat parts and irregular oscillating parts. Time evolution of (V0) 2 at the same 
point is also given in FigfJ(b). One can observe that irregular oscillating parts correspond to 
large (V0) 2 . Further increasing the external forcing, the domain pervades the entire region 
and frequency-locked patterns form. 
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At high u, increasing the amplitude c gives labyrinths or targets. The labyrinths have a 
finite (V0) 2 , which helps the frequency locking now. So, locked labyrinths can arise outside 
the Arnol'd tongue of a single oscillator. But the targets at 3:1 resonance are not frequency 
locked. Although it might be due to the chirality, the origin of unlocked target patterns is 
not clear at this stage and should be studied further. 

In summary, we examined frequency-locking phenomena in spatially extended systems 
using a complex Ginzburg-Landau equation with a resonant term. A resonant tongue struc- 
ture was constructed and studied as a function of //(measuring frequency mismatch) and c 
(amplitude of forcing). Analysis with the phase equation showed that frequency locking is 
enhanced by diffusion at higher forcing frequencies, but suppressed at lower forcing frequen- 
cies. This explains the locus of locked labyrinths and unlocked chaotically bursting domains 
in parameter space with respect to Arnol'd tongues. 

The results are in agreement with recent experimental reports |3|-|5|]. For example, 
labyrinths formed at higher frequencies than tt fronts; 27r/3 fronts at 3:1 resonance move 
due to the chirality. In addition, novel patterns such as chaotically bursting domains and 
target patterns were also observed during the transition to locking. It would be interesting 
to know whether these are confirmed experimentally as well. 

The author wishes to thank H.-T.Moon, S.-O.Jeong and T.-W.Ko for many stimulating 
discussions. This work was supported in part by the project HPC-COSE from the Ministry 
of Science and Technology and in part by BK21 program of the Ministry of Education in 
Korea. 
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frequency of patterns, enhancing locking at low frequencies. 
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FIGURES 

FIG. 1. Phase diagram : frequency deviation v versus critical forcing amplitude c. (a) n = 1, 
(b) n = 2, (c) n = 3. Solid lines are critical amplitude q d for a single oscillator. The shaded regions 
are obtained from the numerical integration of Eq. (|l|) IF: flat, 2Y-.TT fronts, 3F:2-7r/3 fronts, BD: 
chaotically bursting domains, LABY: labyrinthine patterns, TA:targets. See the text for details. 

FIG. 2. Numerical integration results, Plot of Re(^), (a) flat (n = 1, v = -0.4, c = 0.09). (b) 
vr fronts (n = 2, v = -0.1, c = 0.4). (c) 2vr/3 fronts (n = 3, v = -0.1, c = 0.45). (d) labyrinths 
(n = 2, v = 0.1, c = 0.4). (e) target (n = 3, v = — 0.15,c = 0.32). (f) chaotically bursting domains 
(n = 3,i/ = -0.65, c = 0.26). 

FIG. 3. (a) Time series at one point in Fig|2|(c). (b) Parametric plot of real and imaginary 
parts of tj) in a complex plane corresponding to Fig.Q(c). The abscissa is Re(^) and the ordinate 
is Im(^). 

FIG. 4. (a) Chaotic time series at one point in FigJ||(f). (b) Time evolution of (V</>) 2 at the 
same point. 
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This figure "0101003.fig2.jpg" is available in "jpg" format from: 
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